3.2141 \(\int \frac{(a+b x) \sqrt{d+e x}}{(a^2+2 a b x+b^2 x^2)^{5/2}} \, dx\)

Optimal. Leaf size=217 \[ \frac{e^2 \sqrt{d+e x}}{8 b \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^2}-\frac{e^3 (a+b x) \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{d+e x}}{\sqrt{b d-a e}}\right )}{8 b^{3/2} \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^{5/2}}-\frac{e \sqrt{d+e x}}{12 b (a+b x) \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)}-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}} \]

[Out]

-Sqrt[d + e*x]/(3*b*(a^2 + 2*a*b*x + b^2*x^2)^(3/2)) + (e^2*Sqrt[d + e*x])/(8*b*(b*d - a*e)^2*Sqrt[a^2 + 2*a*b
*x + b^2*x^2]) - (e*Sqrt[d + e*x])/(12*b*(b*d - a*e)*(a + b*x)*Sqrt[a^2 + 2*a*b*x + b^2*x^2]) - (e^3*(a + b*x)
*ArcTanh[(Sqrt[b]*Sqrt[d + e*x])/Sqrt[b*d - a*e]])/(8*b^(3/2)*(b*d - a*e)^(5/2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])

________________________________________________________________________________________

Rubi [A]  time = 0.14836, antiderivative size = 217, normalized size of antiderivative = 1., number of steps used = 6, number of rules used = 5, integrand size = 35, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.143, Rules used = {768, 646, 51, 63, 208} \[ \frac{e^2 \sqrt{d+e x}}{8 b \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^2}-\frac{e^3 (a+b x) \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{d+e x}}{\sqrt{b d-a e}}\right )}{8 b^{3/2} \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)^{5/2}}-\frac{e \sqrt{d+e x}}{12 b (a+b x) \sqrt{a^2+2 a b x+b^2 x^2} (b d-a e)}-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}} \]

Antiderivative was successfully verified.

[In]

Int[((a + b*x)*Sqrt[d + e*x])/(a^2 + 2*a*b*x + b^2*x^2)^(5/2),x]

[Out]

-Sqrt[d + e*x]/(3*b*(a^2 + 2*a*b*x + b^2*x^2)^(3/2)) + (e^2*Sqrt[d + e*x])/(8*b*(b*d - a*e)^2*Sqrt[a^2 + 2*a*b
*x + b^2*x^2]) - (e*Sqrt[d + e*x])/(12*b*(b*d - a*e)*(a + b*x)*Sqrt[a^2 + 2*a*b*x + b^2*x^2]) - (e^3*(a + b*x)
*ArcTanh[(Sqrt[b]*Sqrt[d + e*x])/Sqrt[b*d - a*e]])/(8*b^(3/2)*(b*d - a*e)^(5/2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])

Rule 768

Int[((d_.) + (e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Sim
p[(g*(d + e*x)^m*(a + b*x + c*x^2)^(p + 1))/(2*c*(p + 1)), x] - Dist[(e*g*m)/(2*c*(p + 1)), Int[(d + e*x)^(m -
1)*(a + b*x + c*x^2)^(p + 1), x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] && EqQ[2*c*f - b*g, 0] && LtQ[p, -1]
&& GtQ[m, 0]

Rule 646

Int[((d_.) + (e_.)*(x_))^(m_)*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Dist[(a + b*x + c*x^2)^Fra
cPart[p]/(c^IntPart[p]*(b/2 + c*x)^(2*FracPart[p])), Int[(d + e*x)^m*(b/2 + c*x)^(2*p), x], x] /; FreeQ[{a, b,
c, d, e, m, p}, x] && EqQ[b^2 - 4*a*c, 0] &&  !IntegerQ[p] && NeQ[2*c*d - b*e, 0]

Rule 51

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[((a + b*x)^(m + 1)*(c + d*x)^(n + 1
))/((b*c - a*d)*(m + 1)), x] - Dist[(d*(m + n + 2))/((b*c - a*d)*(m + 1)), Int[(a + b*x)^(m + 1)*(c + d*x)^n,
x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && LtQ[m, -1] &&  !(LtQ[n, -1] && (EqQ[a, 0] || (NeQ[
c, 0] && LtQ[m - n, 0] && IntegerQ[n]))) && IntLinearQ[a, b, c, d, m, n, x]

Rule 63

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - (a*d)/b + (d*x^p)/b)^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 208

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-(a/b), 2]*ArcTanh[x/Rt[-(a/b), 2]])/a, x] /; FreeQ[{a,
b}, x] && NegQ[a/b]

Rubi steps

\begin{align*} \int \frac{(a+b x) \sqrt{d+e x}}{\left (a^2+2 a b x+b^2 x^2\right )^{5/2}} \, dx &=-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}}+\frac{e \int \frac{1}{\sqrt{d+e x} \left (a^2+2 a b x+b^2 x^2\right )^{3/2}} \, dx}{6 b}\\ &=-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}}+\frac{\left (b e \left (a b+b^2 x\right )\right ) \int \frac{1}{\left (a b+b^2 x\right )^3 \sqrt{d+e x}} \, dx}{6 \sqrt{a^2+2 a b x+b^2 x^2}}\\ &=-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}}-\frac{e \sqrt{d+e x}}{12 b (b d-a e) (a+b x) \sqrt{a^2+2 a b x+b^2 x^2}}-\frac{\left (e^2 \left (a b+b^2 x\right )\right ) \int \frac{1}{\left (a b+b^2 x\right )^2 \sqrt{d+e x}} \, dx}{8 (b d-a e) \sqrt{a^2+2 a b x+b^2 x^2}}\\ &=-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}}+\frac{e^2 \sqrt{d+e x}}{8 b (b d-a e)^2 \sqrt{a^2+2 a b x+b^2 x^2}}-\frac{e \sqrt{d+e x}}{12 b (b d-a e) (a+b x) \sqrt{a^2+2 a b x+b^2 x^2}}+\frac{\left (e^3 \left (a b+b^2 x\right )\right ) \int \frac{1}{\left (a b+b^2 x\right ) \sqrt{d+e x}} \, dx}{16 b (b d-a e)^2 \sqrt{a^2+2 a b x+b^2 x^2}}\\ &=-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}}+\frac{e^2 \sqrt{d+e x}}{8 b (b d-a e)^2 \sqrt{a^2+2 a b x+b^2 x^2}}-\frac{e \sqrt{d+e x}}{12 b (b d-a e) (a+b x) \sqrt{a^2+2 a b x+b^2 x^2}}+\frac{\left (e^2 \left (a b+b^2 x\right )\right ) \operatorname{Subst}\left (\int \frac{1}{a b-\frac{b^2 d}{e}+\frac{b^2 x^2}{e}} \, dx,x,\sqrt{d+e x}\right )}{8 b (b d-a e)^2 \sqrt{a^2+2 a b x+b^2 x^2}}\\ &=-\frac{\sqrt{d+e x}}{3 b \left (a^2+2 a b x+b^2 x^2\right )^{3/2}}+\frac{e^2 \sqrt{d+e x}}{8 b (b d-a e)^2 \sqrt{a^2+2 a b x+b^2 x^2}}-\frac{e \sqrt{d+e x}}{12 b (b d-a e) (a+b x) \sqrt{a^2+2 a b x+b^2 x^2}}-\frac{e^3 (a+b x) \tanh ^{-1}\left (\frac{\sqrt{b} \sqrt{d+e x}}{\sqrt{b d-a e}}\right )}{8 b^{3/2} (b d-a e)^{5/2} \sqrt{a^2+2 a b x+b^2 x^2}}\\ \end{align*}

Mathematica [C]  time = 0.0243075, size = 68, normalized size = 0.31 \[ \frac{2 e^3 (a+b x) (d+e x)^{3/2} \, _2F_1\left (\frac{3}{2},4;\frac{5}{2};-\frac{b (d+e x)}{a e-b d}\right )}{3 \sqrt{(a+b x)^2} (a e-b d)^4} \]

Antiderivative was successfully verified.

[In]

Integrate[((a + b*x)*Sqrt[d + e*x])/(a^2 + 2*a*b*x + b^2*x^2)^(5/2),x]

[Out]

(2*e^3*(a + b*x)*(d + e*x)^(3/2)*Hypergeometric2F1[3/2, 4, 5/2, -((b*(d + e*x))/(-(b*d) + a*e))])/(3*(-(b*d) +
a*e)^4*Sqrt[(a + b*x)^2])

________________________________________________________________________________________

Maple [B]  time = 0.018, size = 326, normalized size = 1.5 \begin{align*}{\frac{ \left ( bx+a \right ) ^{2}}{24\, \left ( ae-bd \right ) ^{2}b} \left ( 3\,\arctan \left ({\frac{\sqrt{ex+d}b}{\sqrt{ \left ( ae-bd \right ) b}}} \right ){x}^{3}{b}^{3}{e}^{3}+9\,\arctan \left ({\frac{\sqrt{ex+d}b}{\sqrt{ \left ( ae-bd \right ) b}}} \right ){x}^{2}a{b}^{2}{e}^{3}+3\,\sqrt{ \left ( ae-bd \right ) b} \left ( ex+d \right ) ^{5/2}{b}^{2}+9\,\arctan \left ({\frac{\sqrt{ex+d}b}{\sqrt{ \left ( ae-bd \right ) b}}} \right ) x{a}^{2}b{e}^{3}+8\,\sqrt{ \left ( ae-bd \right ) b} \left ( ex+d \right ) ^{3/2}abe-8\,\sqrt{ \left ( ae-bd \right ) b} \left ( ex+d \right ) ^{3/2}{b}^{2}d+3\,\arctan \left ({\frac{\sqrt{ex+d}b}{\sqrt{ \left ( ae-bd \right ) b}}} \right ){a}^{3}{e}^{3}-3\,\sqrt{ \left ( ae-bd \right ) b}\sqrt{ex+d}{a}^{2}{e}^{2}+6\,\sqrt{ \left ( ae-bd \right ) b}\sqrt{ex+d}abde-3\,\sqrt{ \left ( ae-bd \right ) b}\sqrt{ex+d}{b}^{2}{d}^{2} \right ){\frac{1}{\sqrt{ \left ( ae-bd \right ) b}}} \left ( \left ( bx+a \right ) ^{2} \right ) ^{-{\frac{5}{2}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((b*x+a)*(e*x+d)^(1/2)/(b^2*x^2+2*a*b*x+a^2)^(5/2),x)

[Out]

1/24*(b*x+a)^2*(3*arctan((e*x+d)^(1/2)*b/((a*e-b*d)*b)^(1/2))*x^3*b^3*e^3+9*arctan((e*x+d)^(1/2)*b/((a*e-b*d)*
b)^(1/2))*x^2*a*b^2*e^3+3*((a*e-b*d)*b)^(1/2)*(e*x+d)^(5/2)*b^2+9*arctan((e*x+d)^(1/2)*b/((a*e-b*d)*b)^(1/2))*
x*a^2*b*e^3+8*((a*e-b*d)*b)^(1/2)*(e*x+d)^(3/2)*a*b*e-8*((a*e-b*d)*b)^(1/2)*(e*x+d)^(3/2)*b^2*d+3*arctan((e*x+
d)^(1/2)*b/((a*e-b*d)*b)^(1/2))*a^3*e^3-3*((a*e-b*d)*b)^(1/2)*(e*x+d)^(1/2)*a^2*e^2+6*((a*e-b*d)*b)^(1/2)*(e*x
+d)^(1/2)*a*b*d*e-3*((a*e-b*d)*b)^(1/2)*(e*x+d)^(1/2)*b^2*d^2)/((a*e-b*d)*b)^(1/2)/b/(a*e-b*d)^2/((b*x+a)^2)^(
5/2)

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{{\left (b x + a\right )} \sqrt{e x + d}}{{\left (b^{2} x^{2} + 2 \, a b x + a^{2}\right )}^{\frac{5}{2}}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x+a)*(e*x+d)^(1/2)/(b^2*x^2+2*a*b*x+a^2)^(5/2),x, algorithm="maxima")

[Out]

integrate((b*x + a)*sqrt(e*x + d)/(b^2*x^2 + 2*a*b*x + a^2)^(5/2), x)

________________________________________________________________________________________

Fricas [B]  time = 1.09628, size = 1581, normalized size = 7.29 \begin{align*} \left [\frac{3 \,{\left (b^{3} e^{3} x^{3} + 3 \, a b^{2} e^{3} x^{2} + 3 \, a^{2} b e^{3} x + a^{3} e^{3}\right )} \sqrt{b^{2} d - a b e} \log \left (\frac{b e x + 2 \, b d - a e - 2 \, \sqrt{b^{2} d - a b e} \sqrt{e x + d}}{b x + a}\right ) - 2 \,{\left (8 \, b^{4} d^{3} - 22 \, a b^{3} d^{2} e + 17 \, a^{2} b^{2} d e^{2} - 3 \, a^{3} b e^{3} - 3 \,{\left (b^{4} d e^{2} - a b^{3} e^{3}\right )} x^{2} + 2 \,{\left (b^{4} d^{2} e - 5 \, a b^{3} d e^{2} + 4 \, a^{2} b^{2} e^{3}\right )} x\right )} \sqrt{e x + d}}{48 \,{\left (a^{3} b^{5} d^{3} - 3 \, a^{4} b^{4} d^{2} e + 3 \, a^{5} b^{3} d e^{2} - a^{6} b^{2} e^{3} +{\left (b^{8} d^{3} - 3 \, a b^{7} d^{2} e + 3 \, a^{2} b^{6} d e^{2} - a^{3} b^{5} e^{3}\right )} x^{3} + 3 \,{\left (a b^{7} d^{3} - 3 \, a^{2} b^{6} d^{2} e + 3 \, a^{3} b^{5} d e^{2} - a^{4} b^{4} e^{3}\right )} x^{2} + 3 \,{\left (a^{2} b^{6} d^{3} - 3 \, a^{3} b^{5} d^{2} e + 3 \, a^{4} b^{4} d e^{2} - a^{5} b^{3} e^{3}\right )} x\right )}}, \frac{3 \,{\left (b^{3} e^{3} x^{3} + 3 \, a b^{2} e^{3} x^{2} + 3 \, a^{2} b e^{3} x + a^{3} e^{3}\right )} \sqrt{-b^{2} d + a b e} \arctan \left (\frac{\sqrt{-b^{2} d + a b e} \sqrt{e x + d}}{b e x + b d}\right ) -{\left (8 \, b^{4} d^{3} - 22 \, a b^{3} d^{2} e + 17 \, a^{2} b^{2} d e^{2} - 3 \, a^{3} b e^{3} - 3 \,{\left (b^{4} d e^{2} - a b^{3} e^{3}\right )} x^{2} + 2 \,{\left (b^{4} d^{2} e - 5 \, a b^{3} d e^{2} + 4 \, a^{2} b^{2} e^{3}\right )} x\right )} \sqrt{e x + d}}{24 \,{\left (a^{3} b^{5} d^{3} - 3 \, a^{4} b^{4} d^{2} e + 3 \, a^{5} b^{3} d e^{2} - a^{6} b^{2} e^{3} +{\left (b^{8} d^{3} - 3 \, a b^{7} d^{2} e + 3 \, a^{2} b^{6} d e^{2} - a^{3} b^{5} e^{3}\right )} x^{3} + 3 \,{\left (a b^{7} d^{3} - 3 \, a^{2} b^{6} d^{2} e + 3 \, a^{3} b^{5} d e^{2} - a^{4} b^{4} e^{3}\right )} x^{2} + 3 \,{\left (a^{2} b^{6} d^{3} - 3 \, a^{3} b^{5} d^{2} e + 3 \, a^{4} b^{4} d e^{2} - a^{5} b^{3} e^{3}\right )} x\right )}}\right ] \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x+a)*(e*x+d)^(1/2)/(b^2*x^2+2*a*b*x+a^2)^(5/2),x, algorithm="fricas")

[Out]

[1/48*(3*(b^3*e^3*x^3 + 3*a*b^2*e^3*x^2 + 3*a^2*b*e^3*x + a^3*e^3)*sqrt(b^2*d - a*b*e)*log((b*e*x + 2*b*d - a*
e - 2*sqrt(b^2*d - a*b*e)*sqrt(e*x + d))/(b*x + a)) - 2*(8*b^4*d^3 - 22*a*b^3*d^2*e + 17*a^2*b^2*d*e^2 - 3*a^3
*b*e^3 - 3*(b^4*d*e^2 - a*b^3*e^3)*x^2 + 2*(b^4*d^2*e - 5*a*b^3*d*e^2 + 4*a^2*b^2*e^3)*x)*sqrt(e*x + d))/(a^3*
b^5*d^3 - 3*a^4*b^4*d^2*e + 3*a^5*b^3*d*e^2 - a^6*b^2*e^3 + (b^8*d^3 - 3*a*b^7*d^2*e + 3*a^2*b^6*d*e^2 - a^3*b
^5*e^3)*x^3 + 3*(a*b^7*d^3 - 3*a^2*b^6*d^2*e + 3*a^3*b^5*d*e^2 - a^4*b^4*e^3)*x^2 + 3*(a^2*b^6*d^3 - 3*a^3*b^5
*d^2*e + 3*a^4*b^4*d*e^2 - a^5*b^3*e^3)*x), 1/24*(3*(b^3*e^3*x^3 + 3*a*b^2*e^3*x^2 + 3*a^2*b*e^3*x + a^3*e^3)*
sqrt(-b^2*d + a*b*e)*arctan(sqrt(-b^2*d + a*b*e)*sqrt(e*x + d)/(b*e*x + b*d)) - (8*b^4*d^3 - 22*a*b^3*d^2*e +
17*a^2*b^2*d*e^2 - 3*a^3*b*e^3 - 3*(b^4*d*e^2 - a*b^3*e^3)*x^2 + 2*(b^4*d^2*e - 5*a*b^3*d*e^2 + 4*a^2*b^2*e^3)
*x)*sqrt(e*x + d))/(a^3*b^5*d^3 - 3*a^4*b^4*d^2*e + 3*a^5*b^3*d*e^2 - a^6*b^2*e^3 + (b^8*d^3 - 3*a*b^7*d^2*e +
3*a^2*b^6*d*e^2 - a^3*b^5*e^3)*x^3 + 3*(a*b^7*d^3 - 3*a^2*b^6*d^2*e + 3*a^3*b^5*d*e^2 - a^4*b^4*e^3)*x^2 + 3*
(a^2*b^6*d^3 - 3*a^3*b^5*d^2*e + 3*a^4*b^4*d*e^2 - a^5*b^3*e^3)*x)]

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x+a)*(e*x+d)**(1/2)/(b**2*x**2+2*a*b*x+a**2)**(5/2),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [B]  time = 1.2416, size = 463, normalized size = 2.13 \begin{align*} \frac{\arctan \left (\frac{\sqrt{x e + d} b}{\sqrt{-b^{2} d + a b e}}\right ) e^{3}}{8 \,{\left (b^{3} d^{2} \mathrm{sgn}\left ({\left (x e + d\right )} b e - b d e + a e^{2}\right ) - 2 \, a b^{2} d e \mathrm{sgn}\left ({\left (x e + d\right )} b e - b d e + a e^{2}\right ) + a^{2} b e^{2} \mathrm{sgn}\left ({\left (x e + d\right )} b e - b d e + a e^{2}\right )\right )} \sqrt{-b^{2} d + a b e}} + \frac{3 \,{\left (x e + d\right )}^{\frac{5}{2}} b^{2} e^{3} - 8 \,{\left (x e + d\right )}^{\frac{3}{2}} b^{2} d e^{3} - 3 \, \sqrt{x e + d} b^{2} d^{2} e^{3} + 8 \,{\left (x e + d\right )}^{\frac{3}{2}} a b e^{4} + 6 \, \sqrt{x e + d} a b d e^{4} - 3 \, \sqrt{x e + d} a^{2} e^{5}}{24 \,{\left (b^{3} d^{2} \mathrm{sgn}\left ({\left (x e + d\right )} b e - b d e + a e^{2}\right ) - 2 \, a b^{2} d e \mathrm{sgn}\left ({\left (x e + d\right )} b e - b d e + a e^{2}\right ) + a^{2} b e^{2} \mathrm{sgn}\left ({\left (x e + d\right )} b e - b d e + a e^{2}\right )\right )}{\left ({\left (x e + d\right )} b - b d + a e\right )}^{3}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x+a)*(e*x+d)^(1/2)/(b^2*x^2+2*a*b*x+a^2)^(5/2),x, algorithm="giac")

[Out]

1/8*arctan(sqrt(x*e + d)*b/sqrt(-b^2*d + a*b*e))*e^3/((b^3*d^2*sgn((x*e + d)*b*e - b*d*e + a*e^2) - 2*a*b^2*d*
e*sgn((x*e + d)*b*e - b*d*e + a*e^2) + a^2*b*e^2*sgn((x*e + d)*b*e - b*d*e + a*e^2))*sqrt(-b^2*d + a*b*e)) + 1
/24*(3*(x*e + d)^(5/2)*b^2*e^3 - 8*(x*e + d)^(3/2)*b^2*d*e^3 - 3*sqrt(x*e + d)*b^2*d^2*e^3 + 8*(x*e + d)^(3/2)
*a*b*e^4 + 6*sqrt(x*e + d)*a*b*d*e^4 - 3*sqrt(x*e + d)*a^2*e^5)/((b^3*d^2*sgn((x*e + d)*b*e - b*d*e + a*e^2) -
2*a*b^2*d*e*sgn((x*e + d)*b*e - b*d*e + a*e^2) + a^2*b*e^2*sgn((x*e + d)*b*e - b*d*e + a*e^2))*((x*e + d)*b -
b*d + a*e)^3)